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We show that the magnetic Raman response of a spin-liquid is independent of the polarizations 
of the light for triangular symmetries. In contrast, a ground-state that has a broken symmetry 
shows characteristic oscillations when the polarizations are rotated. This would allow to detect 
weak broken symmetries and emergent order-parameters. We focus on the Kagome antiferromagnet 
where no conventional long-range order has been found so far, and present the Raman cross-section 
of a spin-liquid and a valence bond crystal (VBC) using a random phase approximation. 



Interacting spins on the Kagome lattice are particu- 
larly interesting as a possible realization of a spin-liquid 
ground state. The recent discovery of a magnetic oxide, 
ZnCu3(OH)6Cl3 [H, with the geometry of the Kagome 
lattice has triggered renewed interest in this field 0] . AH 
recent experiments in this compound point to the ab- 
sence of long-range order, at least of conventional form, 
and at temperatures well below the superexchange cou- 
pling: a persistent dynamics at the lowest temperatures 
[3|, no Bragg peaks in neutron scattering [5|. This 
is all consistent with the prediction of a non-magnetic 
ground state by exact diagonalizations of the Heisenberg 
model d, S Q- Although other interactions (such as 
that of Dzyaloshinski-Moriya symmetry) may play a role 
in this compound [9], the nature of the ground-state of 
the Heisenberg model is currently under debate. It could 
be a spin-liquid which does not break any symmetry [lo| , 
or a valence-bond-crystal (VBC) that breaks translation 
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symmetry [ll|, |12|, [l3[ [ij, ll5|, ll6|, ll7| . In the latter cases 



however, the unit-cells turn out to be rather large and 
the modulation of the correlations must be weak, other- 
wise a large spin-gap would exist. This makes the use of 
exact spectra to identify the broken symmetry a difficult 
task [l^l , and no definite proof has emerged so far on the 
theoretical side. The question then arises as to whether 
it could be possible to detect such a weak-broken sym- 
metry experimentally. Detecting the resulting emergent 
order-parameters is indeed not easy, for lack of direct 
experimental probes: indirect ways were suggested for 
the chiral order-parameter, for instance, through polar- 
ized neutron experiments [l9| but those do not probe the 
singlet sector of VBC's. 

In this Letter, we argue that the Raman scattering 
of light by spin excitations provides a way to distin- 
guish between a spin-liquid and a broken-symmetry state. 
In particular, it allows one to extract emergent order- 
parameters in systems with no conventional magnetic 
long-range order. Indeed, we have found that the Raman 
cross-section is polarization-independent for a spin-liquid 
governed by a Heisenberg Hamiltonian in triangular ge- 
ometries, whereas it acquires a characteristic polarization 



dependence for a broken-symmetry state, with an ampli- 
tude proportional to the emergent order-parameter. 

Magnetic Raman scattering is a two-photon process: 
the (polarized) electric-field of light forces the exchange 
of two electrons (and hence their spins), thus creating a 
spin excitation in the system. The cross-section of this 
process is given (at zero temperature) by [2^ [2l|, [S^l 

5(w,e,„,eo„0 = I]l(/l^fl|0)|'(5(c^-w/) (1) 
/ 

Hif = ^ {ein ■ ry)(eo„t • ry)Si • Sj (2) 

<^,j> 

^ depends explicitly on the polarizations of the incom- 
ing and outgoing photons, ein,out through the dipolar 
factors e • Tij, where is the bond vector of nearest- 
neighbor sites. The actual polarization dependence re- 
sults from the symmetries: if the ground-state |0) is a 
spin- liquid (that does not break any lattice symmetries), 
one can predict the polarization-dependent terms accord- 
ing to the symmetries of the excited states |/) (of energy 
LOf). For instance, in square geometries, all scattering 
channels depends on polarizations [23| . In triangular 
geometries, however, the cross-section turns out to be 
polarization-independent for a spin-liquid, as is shown 
below using symmetry arguments. Any departure from 
isotropy therefore signals a broken-symmetry state. 

The problem of the energy dependence of the Ra- 
man cross-section ([1]) is also of particular interest: what 
form of singlet response do we expect for a non-magnetic 
ground-state? In a conventional ordered-state with 
(Si) 7^ 0, the response ([T]) is dominated by two-magnon 
excitations at low energy and can be calculated accord- 
ingly [i^. Since this cannot obviously apply to a state 
with (Si) = 0, we have developed the simplest random 
phase approximation (RPA) for the singlet dynamics of 
putative spin-liquid or VBC states. RPA is a rather dras- 
tic approximation that cannot reproduce the large singlet 
sector of the Kagome lattice 0, B] . Nonetheless it gives 
simple well-defined excitations that illustrate the polar- 
ization properties explicited above and which can be fur- 
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ther tested for various lattices. 

First, we discuss general selection rules and the 
polarization-dependent terms of the cross-section using 
symmetry arguments. For this, the Raman operator Q 
is decomposed in irreducible tensors. The point group 
of the Kagome lattice at k = is that of the triangular 
lattice, Csv (there is an additional parity since each site 
is an inversion center, but it does not play any role in 
the following). Cs^ has two one-dimensional irreducible 
representations (IR) Ai and A2 and one two-dimensional 
IR E. The Raman operator does not have a projection 
onto A2 so we have 



Hb = (e„ • eout)OAi +M Oe 



(3) 



where Oai and are irreducible tensors that trans- 
form according to Ai and E respectively. Oai is in fact 
nothing but the Heisenberg Hamiltonian so there is no 
scattering by Ai excited states. This gives a first se- 
lection rule: if the ground state is a spin-liquid (with Ai 
symmetry, the most general case), only the excited states 
belonging to E (they are twice degenerate) are Raman- 
active. The vector M contains the polarization proper- 
ties and is expressed by M = J^ii^in ' u,){eout ■ Ui)u,, 
where the sum runs over the bonds of the unit-cell, 
and iti are the unit-vectors along these bonds ui = 
(1,0), U2 = (-l/2,V3/2), U3 = (-l/2,-^/3/2). M 
can be simply reexpressed with its coordinates M cx 
(cos(6'i„ -I- 0out),sm{9in + 9out)), where 9in and 9out are 
the angles of the polarization vectors of the incom- 
ing and outgoing photons with respect to the x-axis. 
The Raman cross-section ([Ij at zero temperature is re- 
duced to the Fourier transform of (0|i/fl(i)77fl(0)|0) = 
J2a/3M"M'^{0\O%{t)O^E\'^). We use the irreducible de- 
composition of (0|O^(i)O§|0) and that of 

^ ^ . . 

where the ct^'^ are the Pauli matrices. The first term 

belongs to Ai and the last two terms to E, there is no 
projection onto A2. Each term of ([U is multiplied by 
the matrix element (0|O|;(i)O§|0) of the same symme- 
try. This matrix element is zero or not depending on 
the symmetry of the ground state, (i) The ground-state 
is assumed to be a spin-liquid, i.e. it does not break 
any crystal symmetry (the wave- function |0) transforms 
according to, e.g., Ai). In this case the Raman cross- 
section reduces to the Ai terms in the decomposition, 
M2(0|Oij(t) • Oslo). Remarkably, does not depend 
on Gin or Gout- Therefore, the Raman cross-section of a 
spin- liquid is rotationally invariant: 

Spin liquid : S{uj, ej„, eout) = A{uj) 

This is a special property of the triangular symmetry that 
cannot be general: for example, on the square lattice all 



scattering channels are polarization-dependent 2j|. (ii) 
Now suppose that the ground state spontaneously breaks 
a symmetry of the crystal. For instance we could have 
a Neel state or a VBC. In both cases the wave-function 
contains a superposition of degenerate states belonging to 
different IR, |0) = a\Ai)+h\A2)+c\E) . As a consequence, 
the Raman cross-section contains the additional terms 
of Q, _ ^jy2 ^ cos[2(6'i„ + 9out)] and = 
sin[2(0i„ + 9out)] with a prefactor that depends on the 
cross-terms ac*, etc.: 

Broken lattice symmetry : 

) A{uj) + E{uj) cos[2{9in + 9out) + (f>uj] 

The result depends explicitly on 9i„ and 9out through 
the cos term ((f>^ is a phase factor independent of 9i„ 
and 9 out)- The amplitude of the oscillation E{uj) is re- 
lated to the order-parameter of the broken-symmetry 
ground-state and it is weak if the symmetry is weakly 
broken. It appears as a measure of the cross-terms 
of the ground state wave-function. Furthermore if the 
symmetry is broken, we can possibly see individual ex- 
cited states I/) with a transition probability given by 
|(/|M • O£|0)|2 ~ cos2(0„ + 9out + (t>f)- If it is possible 
to resolve individual peaks experimentally then the vari- 
ation with polarizations is strong, otherwise the sum of 
them reduces to the expression above. 

In conclusion, Raman spectroscopy appears as an ap- 
propriate probe to show whether the ground-state breaks 
the crystal symmetries or not: if the response is rotation- 
ally invariant we can conclude that the ground state is a 
spin-liquid [i^ ; if not the ground-state should break the 
symmetry of the crystal. In this case, the amplitude of 
the modulation gives access to the order-parameter. In 
the case of the Kagome system where no ordered moment 
has been found, this may help experimentally to distin- 
guish between a spin-liquid and a VBC for instance. 

This discussion also applies to the triangular lattice. 
Incidentally, the Raman response of the Heisenberg an- 
tiferromagnet on the triangular lattice has been calcu- 
lated recently using exact diagonalizations of small clus- 
ters [2^. It is important to stress that this result of ^ 
cannot be seen as reflecting the Raman response of the 
ordered state in the thermodynamic limit: the oscilla- 
tion has been missed because in a finite-size system the 
ground state belongs to the trivial representation of the 
group. One would need to construct the semi-classical 
Neel state by summing wave-functions of different IR 
with the correct amplitudes. The selection rules and the 
polarization dependence should then reflect the proper- 
ties given above. 

In order to have more precise predictions for the Ra- 
man spectrum, we now present the results of a random 
phase approximation. The approach consists of writing 
down hierarchical equations of motion for the singlet op- 
erator Si{t) ■ Sj{t) (with i and j anywhere on the lat- 
tice) and the four-spin susceptibility associated with ([1]) . 
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The closure of the hierarchy is done in such a way as 
to use two-point correlation functions, (Si • S^) as decou- 
pling parameters instead of the local magnetizations (S^) 
that are assumed to vanish in the system. This involves 
writing down the equation of motion to second-order in 
the time-derivative. The closed equation of motion then 
contains c-numbers that are the Fourier components of 
(Si • Sj). These numbers are unknown and must be de- 
termined self-consistently using sum-rules. This is simi- 
lar to the Kondo-Yamaji decoupling in one-dimensional 
systems 27[, and can be viewed as a random phase ap- 
proximation. For the Kagome lattice with no broken 
symmetry, the self-consistent parameters were previously 
determined 



29| ■ In addition, we can construct a phe- 
nomenological theory of a VBC that breaks the spatial 
symmetries, simply by imposing an ad hoc modulation 
of the (Si • Sj) over the previous parameters. For both 
systems, the equations of motion are solved numerically 
in the reciprocal space and the Raman cross-section ([T]) 
is extracted at k = (the system contains up to 2700 
sites). 
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FIG. 1; (color online). Magnetic Raman cross-section of the 
Kagome lattice calculated within RPA, assuming a spin-liquid 
ground-state. The response does not depend on the orienta- 
tion of the incoming and outgoing photon polarizations, ei„ 
or eout- Each individual delta function is slightly broadened 
with a Lorentzian. 

In Fig. [U we give the result of the RPA calculation 
of S{llj, Bin, Bout) for the Kagome lattice, assuming that 
(Si • Sj) do not break the crystal symmetries (we replace 
the 6 functions in Eq. ^ by Lorentzians). We find that 
the spectrum consists of several peaks, the intensity of 
which decreases strongly when the energy increases. The 
main intensity is in fact quite concentrated in a single 
mode at about 0.8J. The result is fully in agreement 
with the group theory arguments given above: (i) all 
Raman-active modes are twice degenerate, so their wave- 
functions can be labelled with the IR E, (ii) we have 
rotated the photon polarizations ei„ and eout and have 
found that the intensities in Fig. [1] do not change at all. 



By comparing with preliminary exact results of the dimer 
dynamical response of small clusters [s^l , we believe that 
the response in the single peak at w ~ 0.8J can in fact 
be broadened. This would be interesting as an indication 
that the present simple excitation is coupled to other low- 
energy singlet modes and would decay accordingly. Of 
course, the present approach has neglected this effect in 
closing the hierarchy of higher-order Green's functions. 
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FIG. 2: (color online). Magnetic Raman cross-section of the 
Kagome lattice with weak broken-symmetry, within RPA. 
The ground-state is supposed to be the valence bond crys- 
tal depicted in the inset, that breaks the rotation symmetry 
(the correlations are assumed to be stronger on horizontal 
bonds). Different curves correspond to different orientations 
of the polarization vector ei,i(|| Bout)- 

Furthermore, we now choose a VBC ground-state with 
a spontaneous broken symmetry in the two-point corre- 
lation functions. For simple illustration, we impose the 
same weak perturbation of slightly stronger (Si • Sj) on 
all horizontal bonds (30 % stronger in the following cal- 
culation) (see inset of Fig. [J). Since the equation of mo- 
tion within RPA breaks the C^y symmetry, the two-fold 
degeneracy of the ii^-states is lifted and we have pairs 
of peaks. The intensities now display strong orientation 
dependence when ei„ || eout is rotated in the plane, as 
shown by the different curves in Fig. [51 The intensities of 
the two peaks change mainly like cos^ 26 {9in = dout = S) 
for one component and sin^ 29 for the other with dif- 
ferent prefactors as long as the order-parameter is non- 
zero. If the peaks cannot be resolved experimentally, the 
cross-section measures the sum of the two intensities that 
therefore reduces to A{ijJo) + E{ujq) cos 40, in agreement 
with the general argument given above. It should be em- 
phasized that even if the symmetry is weakly broken, the 
sum still displays a weak characteristic oscillation, the 
amplitude of which gives access to the order-parameter 
of the broken symmetry state (sij. We have shown this 
on the simplest example of a VBC, but we believe this 
will remain true for more complicated superstructures 
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with k ^ 0, such as those suggested in the hterature 
ni, [H, [ll Q, [H, [H, [l3| . It is interesting to note that a 



VBC will be presumably accompanied by lattice distor- 
tions of the same symmetry, which can be tested inde- 
pendently by X-rays. 

We conclude with a simple selection rule for Raman 
spectroscopy in triangular geometries: if the ground state 
is a spin-liquid, the Raman response is independent of 
the polarizations of the incoming and outgoing photons. 
However if the ground state has a broken lattice sym- 
metry, it should depend, in most cases [25| . on the po- 
larizations of the light. The dependence is given by 
cos[2{9i„ + 9 out) + with an amplitude that measures 
the strength of the emergent order-parameter. Raman 
spectroscopy can therefore test directly the presence of 
such broken symmetries. This may help in clarifying the 
ground state of the Kagome system ZnCu3(OH)6Cl3 and 
especially to discriminate between a real spin-liquid and a 
valence-bond-crystal. Other particularly interesting can- 
didates are the possible spin-liquids on triangular lat- 
tices: the organic material /«-(BEDT-TTF)2Cu2(CN)3 (it 
is nearly an isotropic triangular lattice) |32| . or NiGa2S4 
[3^ . Here again the issue of emergent broken-symmetry 
is of particular interest. 
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